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ABSTRACT 

We derive a semi-empirical galactic initial mass function (IMF) from observational 
constraints. We assume that the star formation rate in a galaxy can be expressed as the 
product of the IMF, ip(m), which is a smooth function of mass m (in units of Mq ), and 
a time- and space-dependent total rate of star formation per unit area of galactic disk, 
q^T- The mass dependence of the proposed IMF is determined by five parameters: the 
low-mass slope 7, the high-mass slope — r (taken to be the Salpeter value, -1.35), the 
characteristic mass m c h (which is close to the mass m. pea k at which the IMF turns over) , 
and the lower and upper limits on the mass, rri£ (taken to be 0.004) and m u (taken to 
be 120). The star formation rate in terms of number of stars per unit area of galactic 
disk per unit logarithmic mass interval, is proportional to the IMF: 

<*(m) = ^^^rn = S*T^{m) = C^ T m~ T {l - exp [-(m/m ch ) 7+r ] } , 

where AC is the number of stars, mi < m < m u is the range of stellar masses. The 
values of 7 and m c h are derived from two integral constraints: i) the ratio of the number 
density of stars in the range m = 0.1 — 0.6 to that in the range m = 0.6 — 0.8 as inferred 
from the mass distribution of field stars in the local neighborhood, and ii) the ratio of 
the number of stars in the range m = 0.08 — 1 to the number of brown dwarfs in the 
range m = 0.03 — 0.08 in young clusters. The IMF satisfying the above constraints 
is characterized by the parameters 7 = 0.51 and m c h = 0.35 (which corresponds to 
^peak = 0.27). This IMF agrees quite well with the Chabrier (2005) IMF for the entire 
mass range over which we have compared with data, but predicts significantly more 
stars with masses < 0.03 M ; we also compare with other IMFs in current use. We give 
a number of important parameters implied by the IMF, such as the fractional number 
of brown dwarfs and high-mass stars formed at a given time, the average mass of a 
newly-formed star, and the mass of stars formed per high-mass star. 
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1. Introduction 

The Initial Mass Function (IMF) is a fundamental ingredient for modeling any system con- 
taining stars - from clusters to galaxies to the luminous universe. Since the pioneering work by 
Salpeter (1955), the IMF has been derived for a variety of systems, such as clusters of different 
ages, field stars, the galactic bulge, globular clusters, and nearby galaxies. Comprehensive studies 
by Miller & Scalo (1979), Scalo (1986), Kroupa (2001), and Chabrier (2003b, 2005) have provided 
the community with standard IMFs that have allowed the construction of innumerable models with 
a common basis. A comprehensive review of the current state of IMF studies has been given by 
Bastian, Covey & Meyer (2010). Earlier reviews of the theoretical and observational studies of the 
IMF are in Scalo (1998a) and (1998b) respectively. It is impressive that after more half a century, 
the Salpeter work is still the accepted IMF for intermediate and high-mass stars. At low masses the 
situation has evolved much more. The ever increasing sensitivity and resolution of observations, 
together with the theoretical means to infer the stellar masses from the observable quantities, have 
resulted in the discovery of a turn-over in the IMF at ~ 0.5 M and in a progressive reduction of 
the errors in the low-mass IMF. Only few years ago the existence of a flattening or a turnover in 
the IMF at subsolar masses was controversial, but now its existence is well established. However, 
uncertainties in the IMF grow as mass decreases below the hydrogen-burning mass limit (~ 0.08 
Mq ) into the realm of brown dwarfs. 

One of the problems that plagues attempts to determine the IMF is that of unresolved binary 
(or multiple) star systems. One can distinguish three different IMFs: 

• the individual-star IMF, or "true" IMF, ip(m), in which each star that is a member of a 
multiple system is counted separately; 

• the system IMF, ip sys (m), in which each multiple stellar system and each single star is counted 
as one object with its total mass; and 

• the effective system IMF (or effective IMF for short), tfj cS (m), in which unresolved binaries 
are counted as single stars with an effective mass m e ff. Note that Chabrier (2003b) adopted 
a different terminology: he used the term "system IMF" for t/j eS and did not consider what 
we define as ^> sys . 

The relation between the effective mass and the system mass depends on how the mass is deter- 
mined. For stars with known distances (such as those in clusters), the absolute magnitude of the 
system is observed and the effective mass is inferred from a mass-magnitude relation. For field 
stars that do not have independent distance determinations, a color-magnitude relation is used to 
infer the absolute magnitude of the system and then a mass-magnitude relation is used to infer 
the effective mass. Interestingly, Gould, Bahcall, & Flynn (1996) showed that the stellar mass per 
unit area of Galactic disk inferred from the effective IMF agrees reasonably well with that from the 
true IMF. It should be noted that theories of the IMF, such as those of Padoan & Nordlund (2002; 
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see also Padoan et al. 2007) and of Hennebelle & Chabrier (2008, 2009) predict the system IMF, 
whereas observations that cannot resolve close binaries determine an effective IMF. 

The goal of this paper is to infer the individual-star IMF for stars in the disk of the Galaxy. 
In this paper, we shall use the term "star" to refer to star-like objects, whether they burn nuclear 
fuel or not, as well as stellar remnants. Thus, brown dwarfs, white dwarfs, neutron stars and 
stellar-mass black holes (which are sometimes referred to as "collapsed stars") are included in the 
term "star." The terminology for objects below the hydrogen-burning limit has not been fixed yet 
(e.g., see Basri 2003, Chabrier et al. 2005); we use the term "brown dwarfs" for all stars below the 
hydrogen-burning limit. Evolved objects such as white dwarfs do not enter the IMF; the IMF thus 
covers brown dwarfs and main-sequence (i.e., hydrogen-burning) stars. We distinguish stars from 
planets, which form in circumstellar disks and have an excess of heavy elements (Chabrier et al. 
2007, Whitworth et al. 2007). In principle, free-floating planets could be confused observationally 
with brown dwarfs (e.g., Liu & Butler 2002); based on observations of extrasolar planets, such 
free-floating planets are expected to have masses less than 10 Jupiter masses. 

In addressing this problem, we assume that the IMF for stars in galactic disks is both universal 
and simple: 

Universality: Although significant variations in the IMF are observed in the disk of the Galaxy 
(e.g., Scalo 1998a), for the most part these are consistent with random sampling from a universal 
IMF (Elmegreen 1997, 1999; Kroupa 2002; Bastian et al. 2010). Fluctuations about the mean IMF 
can arise in IMF models based on deterministic chaos (Sanchez & Parravano 1999). Evidence for 
evolution of the IMF with redshift is summarized by Elmegreen (2009), but none of this evidence 
deals specifically with disk galaxies. Treu et al. (2010) find that for a sample of early-type galaxies 
a Salpeter IMF provides light-to-mass ratios more consistent with observation than does a Chabrier 
(2003b) IMF, and conclude that massive early-type galaxies cannot have both a universal IMF and 
a universal dark matter halo. Again, this result does not deal with galactic disks. Observational 
evidence that the upper part of the IMF varies with surface brightness is inconclusive: Meurer et 
al. (2009) found that the ratio of Ha to FUV increased with surface brightness, which could be 
explained if massive stars form only in massive clusters (e.g., Pflamm-Altenburg et al. 2009) or if 
the formation of massive stars requires high gas surface densities (Krumholz & McKee 2008). On 
the other hand, detailed observations of outer regions in Arp 78 (Kotulla et al. 2008) and M81 
(Gogarten et al. 2009) indicate that the high-mass IMF is normal in these regions. In any case, 
most star formation occurs in the disk regions where massive stars also form, and we assume that, 
when averaged over large regions of these parts of galactic disks and over long intervals of time, 
the mass distribution of stars is drawn from a universal IMF, ifj(m). This is consistent with the 
conclusion of Bastian et al. (2010), who conclude that there is no clear evidence that the IMF 
varies strongly and systematically as a function of initial conditions after the first few generations 
of stars. 

Simplicity. We further assume that the IMF is a smooth function of stellar mass that is char- 
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acterized by a relatively small number of parameters. There are no known physical processes that 
would make the IMF complex, particularly after averaging over the disparate physical conditions in 
Galactic star-forming regions. The molecular clouds in which stars form are highly turbulent (see 
Vazquez-Semadeni et al. 2000 and McKee & Ostriker 2007 for reviews), and the only scales that 
have been identified observationally in these clouds are the thermal Jeans mass and the mass of the 
entire molecular cloud (Williams, Blitz, & McKee 2000). The close correspondence between the 
mass function of molecular cores and the IMF (reviewed in McKee & Ostriker 2007) suggests that 
this lack of scale should carry over from molecular clouds to stellar systems. There are several ways 
in which complexity might enter into the individual-star IMF: the stellar mass is only a fraction of 
the mass of the core from which it forms, the rest being ejected by protostellar outflows and, for 
massive stars, by photoevaporation; the formation of binaries, etc., is due to fragmentation of the 
core; and finally, the multiplicity of stellar systems evolves with time. In each case, the associated 
physical processes could in principle leave their mark on the IMF. We assume that these effects are 
not large, and in any case they tend to be averaged out by considering the total population of stars 
from many different star- forming regions. 

A significant departure from simplicity would occur if there is more than one distinct star 
formation mechanism. Indeed, several mechanisms have been suggested for the formation of very 
low-mass stars (VLMSs) and brown dwarfs (see the review by Whitworth et al. 2007). Thies and 
Kroupa (2007) have argued that observations of the distribution of the distances of separation of 
binaries supports the existence of a separate mechanism for the formation of brown dwarfs and 
VLMSs. As they point out, there is no reason for the star formation mechanism to depend on 
the hydrogen burning limit, so that there should be a significant mass range in which the two 
mechanisms are both operative. In that case, the total IMF would presumably still be smooth and 
relatively simple. 

Contrary to our assumption of simplicity, some observations do show features in the inferred 
IMF; for example, the mass function of nearby stars measured by Reid et al (2002) has upticks both 
at m ~ 0.1 and at m ~ 1.0 (throughout this paper m is the stellar mass in units of solar masses). 
We assume that these apparent features are due to systematic effects such as uncertainties in the 
mass/magnitude relation and the dependence of this relation on metallicity. 

Based on the assumption that the IMF is simple, we are free to adopt a simple analytic 
expression for it. Our approach in this paper is thus quite different from the classical one of 
inferring the IMF at each mass directly from the data (e.g., Scalo 1986). Furthermore, since our 
form for the IMF is simple, we use integral constraints that are observationally relatively well 
determined to infer the parameters of the IMF. The specific functional form for the IMF that we 
adopt is 

ip(m) oc m~ r {l — exp[— (m/m c h) 7+r ]}. (1) 

which, as we shall show below, is consistent with existing observations. We assume that equation 
([TP applies to the mass range mg < m < m u , and that t/j = otherwise; thus our assumed form for 
the IMF has the form of two smoothly joined, cut-off power laws. Since this function approaches a 
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power law at both low stellar masses (ip oc to 7 ) and at high stellar masses (tp oc m~ r ), we term it 
the Smoothed Two-Power Law (STPL) form for the IMF. It has a total of 5 parameters: the two 
power laws, the characteristic mass, m c h, and the upper and lower limits on the mass, mg and m u . 

Historically, this functional form was first proposed by Paresce & De Marchi (2000) to describe 
the Present Day Mass Function (PDMF) of globular clusters. (They termed this form of the 
mass function the Tapered Power Law, but we prefer to emphasize that there are two power laws 
involved.) De Marchi & Paresce (2001) showed that this form described the PDMFs of galactic 
clusters as well, and that in such clusters the characteristic mass increases with age. De Marchi, 
Paresce & Portegies Zwart (2003) suggested that the PDMF for field stars could be built up by a 
superposition of STPL functions with different values of m c h- At a conference in 2003, Parravano, 
McKee & Hollenbach (2006) proposed the STPL form for the IMF; Hollenbach, Parravano & McKee 
(2005) presented an updated version of this form of the STPL form for the IMF. A similar form 
(but with 7 = 0) has been used by Elmegreen (2006) and Elmegreen et al. (2008). De Marchi, 
Paresce & Portegies Zwart (2010) have given a detailed discussion of fitting STPLs to the mass 
functions of a number of young galactic clusters (ages less than 1 Gyr) and globular clusters over 
the mass range (0.1 — 10)M Q . For the young galactic clusters, they find T ~ 1 and 7 ~ 1.5. The 
characteristic mass increases with age, and they infer an initial value of m c h = 0.15. 

Our approach toward inferring the IMF is complementary to that of De Marchi et al. First, we 
consider a broader mass range, using data from 0.03M Q to over IOOMq; the inclusion of higher mass 
stars leads to a somewhat larger value of T than they found. Second, we rely either on very young 
clusters (to infer the brown dwarf fraction) or observations of local field stars (to fix the shape of 
the IMF for low-mass main sequence stars). The local field stars are the result of the dissolution of 
a very large number of galactic clusters of the type analyzed by De Marchi et al. In our approach, 
we infer the IMF of intermediate-mass stars by interpolation between the IMF of subsolar objects 
and the IMF of massive stars, thereby avoiding distortions of the PDMF due to time variations 
in the star formation rate over the lifetime of intermediate mass stars. Our approach is similar to 
that of Kroupa (2001, 2002), who introduced three and four power-law fits to the IMF; our form 
is simpler in that it has one characteristic mass instead of two or three and is more amenable to 
analytic computation. 

It is customary to fit the low-mass part of the IMF with a log-normal form (e.g., Miller k, 
Scalo 1979; Chabrier 2003b, 2005), which is consistent with the probability distribution functions 
(PDFs) for the density found in simulations of supersonically turbulent, non-self-gravitating molec- 
ular clouds (e.g., Nordlund & Padoan 1999; Ostriker, Stone, &: Gammie 2001). However, as first 
shown by Passot &: Vazquez-Semadeni (1998), an equation of state with an adiabatic index less 
than unity leads to a power-law tail in the density PDF at high densities, which would increase the 
rate of production of brown dwarfs compared to a log- normal PDF (Hennebelle & Chabrier 2009). 
Additional physical effects, such as self-gravity (e.g., Li et al. 2003; Bonnell et al. 2008) and pro- 
tostellar outflows (Li &: Nakamura 2006), also alter the properties of the turbulence. Furthermore, 
the low-mass portion of the mass spectrum of density fluctuations is very difficult to determine 
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computationally due to limitations in resolution. Thus, at the present time, theory offers little 
guidance on the actual form of the low-mass portion of the IMF, and it is an important observa- 
tional goal to determine this. Bastian et al. (2010) point out that if the mid-range of stellar masses 
is described by a log-normal form and if there are power-law tails at both high and low masses, 
then a total of 8 parameters is needed to describe the IMF, which is considerably more complicated 
than the STPL form. As shown in £JH the STPL form of the IMF determined here agrees with 
the Chabrier (2005) form, which is log-normal below 1M , to within 30% for m > 0.03; the IMFs 
differ by more than a factor 2 only for m < 0.017. An observational determination of whether the 
log-normal form extends to the lowest mass stars will be challenging, but it is important in order 
to understand the formation of these stars. 

In this paper (Paper I) we use three constraints to calibrate the STPL form of the IMF: (i) the 
slope of the high-mass IMF, — r; (ii) the ratio Rmk of the number of stars (mainly M dwarfs) in 
the range 0.1 — O.6M to the number of stars (mainly K dwarfs) in the range 0.6 — 0.8Mq; and (iii) 
the ratio i?bd of the number of stars in the mass range 0.08 — 1M© divided by the number of objects 
between 0.03 and 0.08M©. In order to test the accuracy of the STPL form of the IMF, as well as 
those proposed by Kroupa (2002), Padoan & Nordlund (2002), Chabrier (2005), and Hennebelle & 
Chabrier (2008), we compare the IMF shapes with the observations summarized in this paper. 

In Paper II (Parravano, McKee, &; Hollenbach 2010) we carry out tests that depend on the 
normalization of the IMF (i.e., surface density of stars in the Galactic disk) as well as on the shape. 
These tests include the local PDMF, the abundance of white dwarfs and red giants, and the value of 
the surface density of all stars at the solar circle. We then infer the current rate of star formation 
from observations of the ionizing photon luminosity from massive stars in the Galaxy; this rate 
is an average over the last ~ 5 Myr. We compare this current rate with the time-averaged rate 
over the lifetime of the Galaxy (~ 11 Gyr), as determined from observations of the current surface 
density of low-mass stars that have survived this entire period. 

This paper is organized as follows. In §2 we define the basic quantities needed to describe 
the IMF and PDMF and we discuss several different proposed forms for the IMF. In §3, we use 
observational data to determine the shape of the IMF. This shape is compared with observation in 
§4, and finally, in §5 we summarize the results. 

2. Forms for the IMF 

2.1. Basic Definitions 

Let d/V*(m) be the number of main-sequence stars and brown dwarfs with masses between m 
and m + dm. For stars above the hydrogen-burning limit (rribd — 0.075 — Burrows et al 2001), we 
take m to be the Zero Age Main Sequence (ZAMS) mass. We express the PDMF of the main- 
sequence stars and brown dwarfs as dM*/dlnm, which is simply the number of such stars in a 
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logarithmic mass interval. We denote the PDMF per unit area of Galactic disk by 

_ (PK 

<j*(m) = , , . (2) 

It is customary in work on the IMF to express the results in terms of base-10 logarithms, but we 
prefer natural logarithms since they are simpler to use in theoretical calculations. For example, 
Miller & Scalo (1979; hereafter MS79) express the number of main-sequence stars per unit area 
between m and m + dm as (f> ms (rn)d log vn\ in our notation, <^ ms — 0ms 

loge = 0.434^ ms . We 

denote the volume density per unit logarithmic mass interval at the mid-plane by n*o(m). This is 
related to the surface density by £*(m) = 2H(m)n*o(m), where H(m) is the effective scale height 
of stars of mass m. Since the scale height depends on the mass of the star, the PDMF is not 
proportional to the volume density. The corresponding stellar mass densities are £*(m) = mq*(m) 
and p*(m) = mn,(m). 

The differential star formation rate (SFR) by number as a function of mass is given by 
c£/V*(m, t)/d\nm, where N*(m, t) includes only stellar births, not stellar deaths. With our assump- 
tion that the IMF is universal, the differential SFR by number can be expressed as the product of 
the total SFR by number, J\*r(t), and the probability that a star is born with a mass m — i.e., the 
IMF, ip(m): 

dN*{m, t) 



= XrWW, (3) 
a mm 

where 

/ tp(m)dhim = 1; (4) 
J me 

recall that m u is the upper limit on the mass of an individual star and rri£ is the lower limit. Note 
that the assumption of a universal IMF allowed us to write the star formation rate as a separable 
function of mass and time, as did MS79 and Scalo (1986). 

The SFR by number per unit area of Galactic disk is then 

The ratio of the current SFR to the value averaged over the age of the disk, to, is the parameter 
defined by MS79: 

&(*>) - (6) 
where t 

«*t) = ^ [° <Mt) dt. (7) 
*o Jo 

Following MS79, who in turn followed Schmidt (1959), we note that the PDMF is comprised of all 

the stars that are younger than their main-sequence lifetime, r(m), or the age of the Galaxy, and 

is related to the star formation rate by 

rto rto 
c(m) = / c»(t)^M dt = (?,r)^(m) / b(t) dt, (8) 

Jto—T'(m) Jto—r'(m) 
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where r'(m) = min[r(m), to]. Observe that by definition we have 

1 I to 



i J\(t)db = \. (9) 



In our notation, the number of stars ever born per logarithmic mass interval dlnm (Schmidt 1959) 

15 f 

i>{m) / ^ T {t)dt = ip(m){^ T )to = OA3^(\ogm), (10) 

where £(logm) is the notation introduced by MS79. Henceforth, for simplicity we shall generally 
omit the argument t in Jv*(m, t) and related quantities. 



2.2. Previous Forms for the IMF 



Here we describe several forms for the IMF that have been proposed in the literature (cf. 
Kroupa 2002, Bastian et al 2010). 



2.2.1. Salpeter IMF 



The original IMF was proposed by Salpeter (1955) and is the simplest: 

V>sai oc m~ r (me < m < m u ) 

with T ~ 1.35. The mean mass is 



(m) = 



r- 1 



l 



me 
m u 



r-i 



3.5m^, 



(11) 



(12) 



where we have omitted a factor (me/m u ) r <C 1 and where the numerical evaluation is for V = 1.35. 
The mass at which ^g a i reaches a peak (the mode of the IMF) is m pca k = me, which we fix to 
me = 0.21 to fit the ratio Rmk (see §3). This power law could be applied to either individual stars 
or to stellar systems. The primary disadvantage of this form for the IMF is that it does not allow 
for the turnover in the IMF at low masses. 



2.2.2. Scalo IMF (effective) 

Scalo (1986) made a careful assessment of the data available at that time and determined the 
resulting IMF, which he presented in a table in the mass range 0.09 — 63. The mean mass depends 
on the age of the disk, and is in the range (m) = 0.41 — 0.52 for a disk age to = 9 — 12 Gyr. The 
peak of the IMF is m pea k = 0.29. For high-mass stars, this IMF was not corrected for the time 
these stars remain obscured by their parent clouds, nor for the necessity of including the largest 
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OB associations, which are very rare (Parravano et al. 2009), in order to obtain a fair sample of 
massive stars; as a result, this IMF is expected to be deficient in high mass stars. Scalo (1998a) 
proposed a three-power law fit for the IMF in the range 0.1 < m < 100, but we shall not consider 
this form here. 



2.2.3. Kroupa IMF (individual/effective) 

Kroupa (2001, 2002) proposed a four-power law fit to the IMF. This form for the IMF appears 
to be the first in which tp has a power-law behavior at low masses with 7 > 0: 



ipK ot m 11 with 



70= 0.7 ±0.7 (0.01 < m < 0.08) 

71 = -0.3 ±0.5 (0.08 < m < 0.50) 

72 = -1.3 ±0.3 (0.50 < m < 1.00) 

73 = -1.3, -1.7 ±0.3 (1.00 <m). 



(Note that we have modified Kroupa's notation so that ip is in our units of stars per logarithmic 
mass interval. Also note that in the limit of low mass, his 70 is our 7; at high mass his 73 is our 
— T.) Kroupa points out that this IMF is a hybrid, referring to individual stars at low masses 
and stellar systems at higher mass (m^3). Based on the results of Sagar & Richtler (1991), he 
estimates that the individual-star IMF will be somewhat steeper at high masses, with T ~ 1.7 as 
listed above (note, however, that Mafz Apellaniz 2008 finds that binarity has a smaller effect on T). 
The peak of this IMF occurs at m pca k = 0.08. In common with the Chabrier IMF and the STPL 
IMF discussed below, this form of the IMF relates the rate of formation of high-mass stars to that 
of solar-type stars by a simple power law. 

The advantage of this form for the IMF is that it is an approximate analytic representation of 
the data. It is somewhat surprising that the boundary between the lowest two mass ranges occurs 
at the brown-dwarf mass, however, since there is no reason that hydrogen burning, which turns on 
long after the star forms, should affect the mass of the star. Another difficulty with this form is 
that it is a hybrid, referring to single stars at low masses and systems at high masses. 



2.2.4- Chabrier IMFs (individual- star and effective) 



Chabrier (2005) presents both individual-star and system IMFs based primarily on the obser- 
vations of nearby stars by Reid et al. (2002). He finds 



exp 



exp 



(logm - log 0.20) 2 
2 x (0.55) 2 

(logm -log 0.25) 2 
2 x (0.55) 2 



(individual-star) 
(effective) 



(14) 
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for m < 1 and 

V> C oc m- L35±a3 (15) 

for m > 1. These IMFs can be described as a "log-normal/Salpeter" shape. The peak masses 
and the mean masses of MS stars for these IMFs are m pea ^ = 0.20, (m) ms = 0.70 and m pea k = 
0.25, (m) ms = 0.79, respectively, for m u = 120 . 

Earlier, Chabrier (2003b) proposed a log-normal/Salpeter IMF with fripeak = 0.08 for the 
individual-star IMF and 7n pea k = 0.22 for the effective system IMF. As we shall see below, this 
form for the individual-star IMF is not consistent with the data we have analyzed to infer the IMF. 
Henceforth, we shall mean the Chabrier (2005) IMF whenever we refer to the "Chabrier IMF." 



2.2.5. Padoan-Nordlund IMF (theoretical system IMF) 

Padoan & Nordlund (2002) and Padoan et al. (2007) proposed the first IMF based on theory 
that approximately captures the observed behavior of the IMF at both low and high masses. The 
premises underlying this theory are that the mass distribution of gravitationally unstable cores is 
generated by the process of turbulent fragmentation, and that the system IMF is determined by 
this core mass function. The proposed shape of the IMF is (Padoan et al. 2007) 

V'pn (m) oc m~ r / 1 + erf 

where m c h is the Bonnor-Ebert mass in terms of the average density and temperature and a 2 is 
the variance of the gas density PDF. Above 1 Mq , Padoan and Nordlund obtain a power-law 
distribution of dense cores consistent with the high-mass IMF slope. At low masses, the mass 
distribution is determined by the PDF of the gas density, which is assumed to be lognormal. 
Padoan et al. (2007) show that for reasonable values of the physical parameters (turbulent Mach 
number A4q = 10, density n = 10 4 cm -3 , temperature T = 10 K, and slope of the velocity power 
spectrum (3 = 1.9, which leads to T = 1.4, m c h = 1 and a = 1.8), the mass distribution in equation 
(|16p becomes similar to the Chabrier (2003b) effective system IMF, which is similar to the Chabrier 
(2005) effective system IMF. The resultant mode of the IMF and the mean mass of MS stars are 
respectively m pea k = 0.24 and (m) ms = 0.75. In reality, the PN theory addresses the distribution 
of core masses that in principle are larger than the stellar system masses they form. However, if 
the shape of the mass distribution of cores and stellar systems are similar, then equation (I16p with 
a larger m c h gives the core mass distribution. Hennebelle and Chabrier (2008, 2009) assume that 
the core mass is 3 times the stellar system mass. 

While this is an impressive attempt at a theory of the IMF, it has a number of problems 
(McKee & Ostriker 2007, Hennebelle & Chabrier 2008, 2009). The basic idea of their model is 
that fragments are generated by shocks that have a compression that is linear in the shock velocity, 
as expected for shocks with velocities that are large compared to both the sound speed and the 



41n(m/m c h) + o~ 



2 



2V2, 



a 



(16) 
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Alfven velocity; however this assumption is valid only on large scales in large GMCs. They assume 
that the number of fragments scales inversely as the cube of the size, but do not justify this choice 
over several plausible alternatives. Furthermore, their theory neglects the turbulent structure in 
the cores that form massive stars (McKee & Tan 2002, 2003; Hennebelle & Chabrier 2009). 



2.2.6. Hennebelle- Chabrier IMF (theoretical system IMF) 

Hennebelle & Chabrier (2008, 2009) have developed a theory for the IMF based on the Press- 
Schechter (1974) theory for cosmological structure formation. It overcomes each of the problems 
cited above for the Padoan & Nordlund theory. However, in contrast to the Padoan-Nordlund 
theory, it ignores the magnetic field. Like Padoan and Nordlund, they assume that that IMF is 
determined by the mass function of the molecular cores out of which the stars form. Since they 
do not follow fragmentation of the cores, their IMF is a system IMF, which they take to be the 
same as the effective IMF. They assume that the stellar system mass is 1/3 of the core mass based 
on observations suggesting that the core masses are several times greater than the resulting stellar 
masses (e.g., Alves, Lombardi & Lada 2007). For an isothermal equation of state, their core mass 
function is (Hennebelle Sz Chabrier 2008) 



(mi? 3 )V2 
where 



r_N&w_£\ (17) 



1 + (1 - r])MlR 2r > 
1 + (2n + l)M 2 R 2l > 

fh= =R(l + MlR 2v ), (18) 

m c h is the Jeans mass, R is the radius of the clump in units of the Jeans length, a 2 is the variance 
of the gas density fluctuations, is the characteristic Mach number at the Jeans scale, and r\ 
is the exponent of the linewidth-size relation. They take n to be in the range 0.4 — 0.45 based on 
the simulations of Kritsuk et al (2007). At moderately high masses (m> m c h), m/R 3 oc m~ os to 
m~ - 7 and the mass dependence due to the expln 2 factor is weak; if that is ignored, 

V>hc oc m~^+t)/{2r,+i) = ^-r HC (m ^ 1), (19) 

which gives Thc — 1-3 for the value of r\ they adopt. At very high masses (m S> 1), the exponential 
term again takes over the mass dependence and -0HC drops off more steeply with mass, reverting 
to a log- normal type distribution. At low masses (rh < 1), m/R? oc m~ 2 and the exponential 
term rises rapidly with increasing mass, producing a log normal form again. Their results for 
a non- isothermal equation of state are considerably more complicated, but yield results that are 
qualitatively consistent with the isothermal theory (Hennebelle & Chabrier 2009). 

Both the Padoan-Nordlund IMF and the Hennebelle-Chabrier IMF apply to a particular star- 
forming cloud. In order to determine an average IMF that can be compared with observations 
of stars from different clouds, it is necessary to average the IMF over a distribution of cloud 
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temperatures, densities and Mach numbers. We assume that it is possible to represent this average 
by a suitable choice of m c h, M.* and a. For the parameter values typically used by Hennebelle and 
Chabrier (77 = 0.4, .M* = v2 and a = 1.9), the mode of the distribution matches the Chabrier 
effective value m pea k = 0.25 for m c h = 1.56. In this case, the mean mass of MS stars is {m) ms = 0.88. 



2.3. Smoothed Two-Power Law IMF 

Existing observations of the IMF are consistent with, but do not require, power-law behavior at 
low masses as well as at high masses. Figure [U shows the Hillenbrand (2004) compilation of several 
studies providing the IMF slope as function of the mass covered by the study. As she pointed 
out, the relatively small scatter at low masses is quite surprising given the variety of observational 
techniques and the variety of transformations between the observations and stellar masses that are 
used. Above masses of several M , the scatter is larger, but the mean slope does not appear to 
change with mean mass and is centered on the Salpeter (1955) slope of -1.35. The Salpeter slope 
holds down to at least 1 M and believably down to 0.6 M . 

To capture a power-law behavior at low masses (ijj oc to 7 ) and high masses (tp oc m~ T ), we 
adopt a simple analytic form, the Smoothed Two-Power Law (STPL) IMF, as discussed in §1. This 
form for the IMF has a shape that is determined by three parameters plus a lower and upper mass 
cutoff: 

-0STPl("i<) = Cm~ r {l - exp [-(m/m ch ) 7+r ] } (m e < m < m u ). (20) 

Each of the five parameters (mg, m u , m c h, 7, and T) has a direct physical significance; the 
normalization factor C is determined in terms of these parameters (see below). We assume that 
star formation is suppressed below a mass ni£ (possibly due to opacity effects in fragmentation — 
Low & Lynden-Bell 1976) and above a mass m u . The parameter m. c h approximately determines 
the position of the IMF maximum. Finally, 7 and — r are the low-mass and high-mass slopes of the 
IMF, respectively. The advantages of this form for the IMF are that it captures the observed power- 
law behavior at low and high masses and that it has a simple analytic form. The disadvantage is 
that, unlike a log-normal, for example, there is no parameter that independently determines the 
width of the IMF; instead, the full-width half maximum of the IMF is determined by the two power 
laws. Nonetheless, as we shall see, it provides a good fit to the available data. 

Since the STPL IMF has a simple, analytic form, it is possible to directly calculate a number 
of quantities of interest: The normalization factor C, which ensures that f™™ V'STPL(m) cilnm = 1, 
is 

rrru. 



1 

c 



I m~ r {l — exp [— (m/m c h) 7+r ] } dhxm, (21) 

J mo 



1 m e 
1 



rm ch 



G 



7 \ r / me 



(22) 



.7 + 17 7 V«Mi, 

where G(x) is the Gamma function and we have assumed that mi <C m c h *C m u . The mean mass 
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Fig. 1. — Hillenbrand (2004) compilation of several studies providing the IMF slope as function of 
the mass covered by the study. The horizontal thin line through each data point is the range of 
masses over which the single slope is claimed to apply. The curves superimposed on the Hillenbrand 
(2004) plot correspond to the mass dependence of the slope of two IMF functions: The continuous 
curve is the t/'stpl (eq. [20] with T = 1.35, 7 = 0.51 and m c h = 0.35), and the dashed curve gives 
the slope of ipc (eqs. fl~4"l[T5|h 
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while the peak mass (accurate to 4% for T > 1 and r/7 > 1.2) is, 
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The fraction of stars that are born as brown dwarfs (i.e., with masses less than TObd) is 



V m ch 
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(25) 



which is accurate to better than 2%, and the mass fraction of stars born as brown dwarfs is 
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mip dlnm, 
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(26) 



(27) 



which is accurate to better than 



2.3.1. De Marchi-Paresce-Portegies Zwart IMF (effective IMF) 

As discussed in the Introduction, Paresce & De Marchi (2000), De Marchi & Paresce (2001), 
and De Marchi, Paresce & Portegies Zwart (2003, 2010, hereafter MPPZ) have applied the STPL 
form of the IMF to both young galactic clusters and globular clusters. They make no attempt to 
correct for binarity, so theirs is an effective system IMF. They found that the characteristic mass of 
a star in a cluster increases with time due to dynamical evolution. Over a mass range 0.1 < m < 10, 
De Marchi et al (2010) infer T = 0.97 ± 0.17, 7 = 1.54 ± 0.64, and, at birth, m ch = 0.15. 

The shape of this IMF and the others considered in this paper are shown in Figure [2] for the 
parameter values given in Table 2. 
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Fig. 2. — Different forms of the IMF normalized to V'STPL- a) Individual-star IMFs ^Sal> ^Cca 
and ipK normalized to VsTPL,(r = 1.35,7 = 0.51, = 0.35). ipc 3 refers to the individual-star 
IMF proposed by Chabrier (2003a, 2003b), which has the same form than tpc but with m pe ak = 0.08 
and a = 0.69. b) Effective system IMFs Vsai> ipc, ipK, V'mppz and tpS86 normalized to V'STPL (r = 
1.35,7 = 0.57, m ch = 0.42). Vsse is the IMF in Scab (1986) Table 7 for a disk age of 12 Gyr and 
&o = 1- c) Theoretical system IMFs Vpn and ipnc normalized to V'stpl(F = 1.35,7 = 0.57, m c h = 
0.42). For all IMFs f ipdlnm = 1 in the range 0.004 < m < 120, except for ^Sai and ipsse for which 
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3. Observational Determination of IMF Parameters (Individual-Star) 

3.1. The High-Mass Slope T 

The STPL form for the IMF has a single value of T for intermediate and high-mass stars. The 
value of r for clusters and associations shows considerably scatter but the average value of T is 
close to the Salpeter value. For clusters, Scalo (1998a) concluded that T for 1-10 M is 1.5 to 1.8, 
but probably flatter (~ 1.2) at larger masses. A fit to the high-mass end of Scalo's (1986) IMF 
gives 1.5 < r < 1.7. However, one obtains 1.3 < T < 1.5 when Scalo's (1986) IMF is corrected by 
considering that (i) short-lived stars spend part of their lifetime obscured by their parent clouds 
and (ii) in order to obtain a fair sample of massive stars, it is necessary to include the largest OB 
associations, which are very rare and therefore often omitted in Galactic surveys (Parravano et al. 
2003, 2009). Additionally, Elmegreen and Scalo (2006) pointed out that if the SFR were decreasing 
because we are moving farther from the density wave shock, but one assumed that it was constant 
in determining the IMF from the PDMF, then one would get a value that is too steep. This effect 
adds to the above two corrections in the same direction. 

When the various estimates of the IMF slope derived from star counts of clusters and associ- 
ations are plotted as function of the average mass over which a common slope is claimed to apply 
(Scalo 1998a, Kroupa 2001, Hillenbrand 2004; see Fig. 1), the data show considerable scatter for 
m > 1, but the average value of T is close to the Salpeter value of 1.35. Elmegreen (2006) reviewed 
observations of galaxy-wide IMFs and concluded that, with rare exceptions, the high-mass slope 
averaged over entire galaxies is consistent with that observed in star clusters, which is about the 
Salpeter value. More recently, Bastian et al (2010) concluded that the available data are indeed 
consistent with V ~ 1.35. De Marchi et al. (2010) find T ~ 1 for young clusters in their sample, 
but it must be borne in mind that they are focusing on stars with m < 10. 

How does binarity affect the high-mass slope? Mafz Apellaniz (2008) has summarized the 
evidence that most massive stars are in binaries. However, through numerical experiments, he 
shows that this does not have a significant effect on the high-mass slope of the IMF. For our 
standard model we therefore adopt T = 1.35, but we explore the consequences of changing T to 1.2 
and 1.5. 



3.2. Upper Mass Limit m u 

The existence of an upper limit to the mass function has been controversial; for example, 
Massey & Hunter (1998) suggested that the apparent cutoff is due to the scarcity of very massive 
stars. This issue was addressed in a general way by MW97 and Parravano et al. (2003), who showed 
that the cutoff is physically significant for various samples of stars in the Galaxy. Elmegreen (2000) 
argued that the data support the existence of an upper mass limit and proposed a model that shows 
a high-mass turndown in the IMF. Oey & Clarke (2005) showed that observations of a number of 
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star clusters indicate that m u < 200. In a detailed study of the Arches Cluster, Figer (2005) 
concluded that m u < 150. These studies do not take multiplicity into account; for example, a star 
that was thought to violate this upper mass limit, Pismis 24-1, with an inferred mass exceeding 
200M Q , has been resolved into three stars with masses in the range 92 — 97M (Mafz Apellaniz 
et al. 2007). The highest dynamically determined stellar masses are those in a Wolf-Rayet binary, 
with masses of 83 ± 5M Q and 82 ± 5M Q (Bonanos et al. 2004); note that these are present-day 
masses, not initial ones. Maiz Apellaniz (2008) summarizes the evidence that the upper mass limit 
is m u ~ 120. 

For our standard model we therefore adopt m u = 120 for individual stars, but we explore the 
consequences of adopting values from m u = 80 — 1500 In reality, the decrease in the IMF near 
the upper mass limit is likely to be smooth rather than abrupt, so that our model of an IMF with 
a sharp cutoff at m u is approximate. The global effects of massive stars are often inferred from 
their ionizing luminosity, but the existence of a few stars with masses exceeding m u in a cluster 
or galaxy with a very large number of massive stars does not have a significant effect on the total 
ionizing luminosity. As stated above, massive star formation might be suppressed in the outer parts 
of galaxies, but most of the star formation occurs in the disk regions where massive stars also form. 

3.3. Lower Mass Limit mi 

Theoretically, opacity effects should limit fragmentation into objects smaller than some critical 
value (Low & Lynden-Bell 1976); including the effects of He, this is about O.OO4M (Whitworth 
et al. 2007). Caballero et al. (2007) have found that the mass distribution of brown dwarfs in one 
cluster extends smoothly down to O.OO6M0, the lowest limit observed to date, which is consistent 
with this theoretical limit. We adopt mi = 0.004. Because the IMF is rising with mass at very low 
masses, the exact value of this lower limit does not appreciably affect our estimates of the number 
or mass fraction of brown dwarfs. 

3.4. The Shape of the Low-Mass Main-Sequence IMF 

To determine the characteristic mass, m c h, and the low-mass power law, 7, we first consider 
the shape of the low-mass, main sequence part of the IMF. We focus on stars that have sufficiently 
low masses (m < 0.8) that their lifetimes exceed the age of the disk; thus, the IMF is the same as 



A very recent paper by Crowther et al. (2010), utilizing observations of the young star clusters NGC3603 and 
R136, suggests that m u may be as large as 300. As we will show, large values of m u have very little effect on the 
conclusions of this paper, and we maintain our range 80< m u < 150 here. 
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the PDMF for these stars. We use one parameter to describe this part of the observed IMF, 

K{m = 0.1-0.6) 
RMK = N*(m = 0.6-0.8)' (28) 

which is very roughly the ratio of the number of M stars to K stars. 

The low-mass, main-sequence IMF can be inferred using local field stars, distant field stars, or 
cluster stars. Each data source has its advantages. We shall use field stars in the local neighborhood, 
which have the advantages that many of the binaries are resolved and the distances are directly 
measured, so that the absolute magnitude is known. The dominant uncertainty in this case is 
the relation between the mass of a star and its absolute magnitude in whatever photometric band 
used for the survey. Samples of distant field stars can be far larger (e.g., Zheng et al. 2001; 
Covey et al. 2008; Bochanski et al. 2010), but have additional disadvantages: First, a metallicity- 
dependent color-magnitude relation must be used to estimate absolute magnitudes and distances, 
and second, one must statistically correct for binaries to estimate the individual-star mass function. 
Clusters have the advantage that all the stars are at the same distance. However, they may have 
systematic deviations from the average IMF due to the circumstances of their formation and again, 
a correction for binaries must be made (nonetheless, we shall use observations of clusters to estimate 
the contribution of brown dwarfs to the IMF-see §3.51 below) . Local Galactic clusters are generally 
not large and so are subject to statistical uncertainties; globular clusters provide large samples of 
stars, but are subject to evolutionary effects (e.g., De Marchi et al. 2010). 

We use the local sample of field stars created by Reid et al. (2002, hereafter RGH02) based on 
data from the Hipparcos catalog and the Palomar /Michigan State University (PMSU) survey. The 
PMSU volume-complete sample (16 > My > 8) includes 558 main-sequence stars in 448 systems 
with distance limits ranging from 5 pc {My = 15.5) to 22 pc {My = 8.5). The Hipparcos 25 pc 
sample {My < 8) includes 538 single stars, 204 binaries, 22 triples and four quadruple systems, 
corresponding to 1028 stars in 768 systems. The luminosity function created from these data 
includes a correction for undetected binaries in the Hipparcos sample: RGH02 find a multiplicity of 
only 30% in that sample, about half that Duquennoy &: Mayor (1991) inferred for the multiplicity 
of solar-type dwarf stars, and they therefore doubled the contribution of the measured secondaries. 
This correction has a maximum effect of 18% in a single magnitude bin, and the typical effect is 
less than 12%. 

We opted not to use the recent very large sample of distant field stars by Bochanski et al. 
(2010) for the reasons just cited as well as the following. They analyzed about 15 million low- 
mass stars (0.1 < m < 0.8) from the Sloan Digital Sky Survey and estimated the mass function 
of field low-mass dwarfs. However, we note that close to the upper mass limit, the adopted mass- 
Mj relation (Delfosse et al. 2000; hereafter D00) needs to be extrapolated for Mj < 5.5 (i.e. 
m e g > 0.72). In addition, the binary correction used to estimate the single-star from the system 
mass function does not account for the secondary stars in systems with effective mass above 0.8 
Mq . For these reasons the individual-star mass function estimated by Bochanski et al. (2010) is 
likely to be an underestimate above 0.7Mq. This is probably the reason why the lognormal fit to 
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this mass function has about the same peak mass, m pea k ~ 0.18 ± 0.02, as Chabrier (2005, see eq. 
[Tl"]) but a smaller dispersion, a = 0.34 ± 0.05 instead of 0.55. Compared with the Pleiades mass 
function (Moraux et al. 2004), the two mass bins above m ~ 0.6 in the Bochanski distribution (see 
Fig. 27 in Bochanski et al. 2010) are about a factor 0.5 below the Pleiades distribution. Therefore, 
it is not possible infer an accurate value for Rmk from the results of Bochanski et al. (2010). In 
§5.1 we compare the effective mass function of distant field stars to our estimation of the effective 
IMF and discuss some possible sources of discrepancy. However, we emphasize again here that the 
main goal in this paper is to constrain the individual-star IMF. 

To determine Rmk an d its uncertainty, we use the luminosity function found by RGH02 for 
a local sample of field stars together with the two mass-luminosity relations they used to derive 
the mass distribution. One of the m{My) relations is based on Delfosse et al. (2000) for My > 10 
and on their own empirical m(Mv) relation for My < 10 [hereafter denoted as the RGH m(My) 
relation]; for this relation, masses 0.1, 0.6 and 0.8 Mq correspond to My = 16.4 , 8.84 and 7.0 
magnitudes, respectively. The other mass-luminosity relation used by RGH02 is the Kroupa, Tout, 
&; Gilmore (1993) semi-empirical mass-luminosity relation [hereafter, the KTG m{My) relation]; 
for this relation, masses 0.1, 0.6 and 0.8 M Q correspond to My = 16.13, 8.7 and 6.24 magnitudes, 
respectively. These two m(My) relations give about the same number of stars in the mass range 
0.1-0.6 Mq , but the KTG relation includes more stars than the RGH one in the mass range 0.6-0.8 
M . With the KTG m{My) relation, the RGH02 luminosity function implies Rmk — 8, whereas 
for the RGH m(My) relation, Rmk — H- We assume that Rmk is between these two values, so 
that 

Rmk = 9.5 ± 1.5 (29) 

for the IMF of individual stars. 

3.5. Ratio of Low-Mass Main- Sequence Stars to Brown Dwarfs 

Uncertainties in the IMF grow as stellar masses decrease below the hydrogen-burning mass limit 
(~ 0.08 Mq ) into the realm of brown dwarfs. The derivation of the mass function of field brown 
dwarfs from a survey in a given band involves the use of theoretical cooling curves parameterized 
by mass and age (Burrows et al. 2001; Baraffe et al. 2003), and a star formation history (Chabrier 
2002). Furthermore, due to the intrinsic low luminosity of evolved brown dwarfs, the numbers of 
objects in these surveys are relatively small. On the other hand, very young embedded clusters 
are particularly valuable because young brown dwarfs have relatively high luminosities and can be 
readily detected using infrared photometry. However, mass functions derived by placing the young, 
low-mass stars on the theoretical H-R diagram using spectroscopic and photometric observations 
are sensitive to the theoretical pre-main sequence (PMS) model used to infer the masses and ages 
of the stars. Furthermore, it is often assumed that all the stars in a cluster have the same age, 
whereas it is possible that the star formation extends over several dynamical times (Tan, Krumholz, 
& McKee 2006). Despite these uncertainties, nearby young clusters provide the best opportunity to 
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determine the IMF in the low-mass and brown-dwarf regime and test whether it might be universal. 

We note that some theories suggest that the low-mass IMF might not be universal: In the 
Padoan & Nordlund (2004) theory of brown-dwarf formation, for example, brown dwarfs form in 
rare, high-density peaks, and therefore their abundance could be exponentially sensitive to the 
ambient conditions such as the Mach number of the supersonic turbulence (McKee & Ostriker 
2007). It is thus of great interest to determine the IMF at very low masses. 

Observations have established that the IMF peaks at subsolar masses and declines into the 
brown-dwarf regime. For example, Luhman et al. (2000) compared the IMFs derived from the 
H-R diagram for the Trapezium cluster by using two different PMS models. With the D'Antona & 
Mazzitelli (1997) PMS model, the IMF has a turnover at ~ 0.3 M , whereas with the Baraffe et 
al. (1998) model, the turnover occurs at ~ 0.6 M . Muench, Lada, & Lada (2000) and Muench 
et al. (2002) derived the Trapezium IMF from the infrared (K-band) luminosity function and 
found a turnover around 0.15 M Q . Despite the quantitative differences between these IMFs, they 
all show that the mass distribution dN*/dlnm has a broad maximum at subsolar masses, with 
a turnover between 0.1 and 0.6 M & , and then declines with decreasing mass in the brown-dwarf 
range. Bastian et al (2010) conclude that the existing data show that the IMF peaks in the mass 
range (0.2 — O.3)M and then declines with 7 > 0.5. 

The fraction of stars that are brown dwarfs sets a strong constraint on the IMF that is relatively 
insensitive to the effects of uncertainties in the PMS tracks. This fraction can be determined 
from observations of nearby young star clusters, which provide samples that extend well below 
the hydrogen-burning limit. Since these data do not take binarity into account, they provide a 
constraint on the effective IMF. We follow Andersen et al. (2006, 2008) in characterizing the low- 
mass IMF function by the ratio of the number of subsolar main-sequence stars to that of readily 
observable (m > 0.03) brown dwarfs in nearby clusters, 

s _ Af*(m = 0.08-1) 

bd "AA,(m = 0.03 - 0.08)' W 

Andersen et al. (2006, 2008) approximated the transition from the main sequence to brown dwarfs 
as occurring at = 0.08. We adopt this value for evaluating R^, but for all other properties 
of the population of brown dwarfs, we use the best theoretical value, mbd = 0.075 (Burrows et al. 
2001). Note that this ratio is not influenced by possible variations of the very low cutoff mass of the 
IMF because m = 0.03 is significantly above the expected values of this cutoff (~ 0.004 — O.OO6M0; 
Whitworth et al. 2007, Caballero et al. 2007). 



3.5.1. Universality of the Brown-Dwarf IMF 

In order to use the ratio R^ d to determine the IMF, it must have a universal value. Andersen 
et al. (2006) concluded that the brown-dwarf ratios they analyzed, including those from field-star 
PDMFs, are all consistent with a single value to within the errors. Andersen et al. (2008) showed 
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that the values of R^ d for seven clusters (Taurus, ONC, Mon R2, Chamaeleon, Pleiades, NGC 2024 
and IC 348) are statistically consistent with the Chabrier (2005) lognormal IMF, a result that is 
consistent with a universal IMF in the brown-dwarf regime. Here we strengthen this conclusion 
by adding more cluster data and by focusing on the universality of R^, not on a particular IMF 
shape. 

In addition to the seven clusters considered by Andersen et al. (2008), we include data from 
recent studies of NGC 1333 (Scholz et al. 2009) and NGC 6611 (Oliveira et al. 2009). Table 1 
gives the brown-dwarf ratios, -Rfo bs , and the number of objects, A/i j0 bs (0.03 < m < 1), observed 
in each cluster, labeled Note that A/^obs is the number of detected objects in the surveyed 
area, not the actual number of objects in cluster i. The estimated errors in -R^Q^g agree with the 
expected root J\f statistical fluctuations in the number of brown dwarfs (A/bd) hi t ne mass range 
0.03 < m < 0.08 and in the number of subsolar stars in the mass range 0.08 < m < 1. 
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Table 1. Ratio of Subsolar Stars to Brown Dwarfs in Young Clusters 



Cluster 
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-^i,obs 


ocff 
1 l i,obs 
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185 


O O + 0.8 

^•^-0.7 


Pleiades 
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200 
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NGC 2024 
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IC 348 
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NGC 1333 


8 
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NGC 6611 
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160 
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The value of the brown-dwarf ratio corresponding to the total observed sample is 



R\ 



.off 

'bd,tot 




= 4.6 



+0.55 
'-0.47' 



(31) 



where a root J\f statistical error was assumed in the numerator and the denominator of equation 
(|30l) . We note that this result is consistent with the Chabrier (2005) effective IMF, which has 

neS _ e 

-"bd — °- 

To determine if the nine values of Rf^ s in Table 1 are consistent with a single universal value, 
-^-bduniv' ana o hi particular, if they are consistent with the sing le value Rf dtoV we compare the 
observed sample {-Rbdi°bs i n Table 1 to 1000 synthetic samples {-Rbdisyn, each one consisting of 
nine synthetic clusters that are generated by imposing statistical fluctuations on the observed ones. 
Let Ri t k be the value of i?bd for cluster i in the synthetic sample {R\^}k, where 1 < i < 9 and 
1 < k < 1000. The values of Ri^ are generated by Poisson fluctuations on the number of subsolar 
stars and on the number of brown dwarfs observed in cluster i. That is, 



where x is the assumed universal value of R^, and J- in) is an integer number drawn from a Poisson 
distribution of mean equal to n. 

The value Rf dxmiv , or x, is varied from 3 to 6.7 in steps of 0.0025. For each value of x 
the following steps are performed: (1) 1000 synthetic samples {R\^}k are generated, where each 
synthetic sample k consists of nine it^d va l ues generated with equation (|32p . (2) The significance 
level P({R^}) of the Kolmogorov-Smirnov test is calculated for the observed sample {R^}oba 
in Table 1 and for each of the 1000 synthetic samples, {R\^}k- A cumulative distribution of a 
Gaussian centered at x and with standard deviation a = 0.3 is used in the Kolmogorov-Smirnov 
test. (3) The values of P ({R\^}k) for the 1000 synthetic samples are ranked in descending order 
(4) The position of the significance level for the observed sample -P({-Rbd}obs) i n the above ranking 
list is obtained. 

Figure [3] shows how the observed sample ranks among the Kolmogorov-Smirnov tests of the 
1000 synthetic samples (the rank is 1 for the best significance level) for each value of x. The 
noise-like fluctuations reflect the random changes in the 1000 synthetic samples. The asymmetry 
of rank (R^ un ; v ) is due to the fact that larger Poisson fluctuations occur in the number of brown 
dwarfs (for R^ > 1), and these fluctuations increase with increasing R^ univ - ^ n an y case > there 
is a pronounced minimum close to ^dtof suggesting that ^bduniv — ^bdtot ^ s highly probable, 
especially because the rank is in the top 20%. This method does not allow us to quantify the error 
in the value of -Rbduniv- However, as shown in Figure O we can quote for example the range of 
-^bduniv ^ or w hich the observed sample ranks in the top 66%; that is, 3.95 < ^duniv < 5.1. Our 
results are fully consistent with those of Andersen et al. (2008): They found that 63% of their 
simulations had a probability higher than what is found for a Chabrier IMF, which has R^ = 5, 



Ri,k 



^V^obs x/{l + x)] 

rWi,<*»/(i + x)] 



(32) 
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while the results in Figure [3] show that the rank of the observed sample for it^d = 5 is ~ 650/1000. 
We conclude that the mass distributions in the observed sample are consistent with a universal 
value of the brown-dwarf fraction 



3.5.2. The Brown-Dwarf Ratio RJj°$ for the Individual- Star IMF 

To infer the individual-star IMF from the constraints Rmk, R^d and T, it is necessary to 
correct the ratio i?bd for unresolved binaries. The ratio -RJJJk = 9.5 ± 1.5 obtained in ^3,41 was 
estimated from the mass function of nearby individual stars, whereas the ratio R^ univ = 4.6+065 
estimated in §3.51 was obtained from observations of young clusters, in which most binary systems 
are unresolved. In Appendix A we make a series of assumptions to account for binaries and then 
allow for the difference between the effective and the actual masses of stars to find (see eq. IA12|) 
that i?jf d d ~ lARf d /(1.2 + gi?bd)- Therefore, the effective brown-dwarf ratio 3.95 < Rf d < 5.1 
corresponds to 3.37 < R^ < 4.04, so that 

K£ = 3-76^:1. (33) 

(Note that in obtaining the range of -Rj^f, we ignored the errors in the parameters in equation lA121 ) 
Although Chabrier's (2005) fit to the IMF focused on stars with m > 0.1, our estimate for -R^dumv 
agrees with the value for the Chabrier effective IMF, Rf~^ ~ 5, and our estimate for i? bd univ agrees 
with the value for the Chabrier individual-star IMF, — 4. Chabrier's values for i?bd he near 
the upper end of the values we infer from the data. 



3.6. Results: Observational Determination of m c h and 7 

To determine the best-fit values of of m c ^ and 7, we calculate level curves for i? b ^ and -RJ^K 
in this parameter space under the assumption that V = 1.35; the values of mg and m u do not 
have a significant effect. Figure Q] shows where ^stpl fulfills the constraints R 1 ^ = 3.761q39 and 
i?MK = 9.5 ± 1.5. For the standard STPL form of the individual-star IMF, we therefore adopt the 
values 7 = 0.51 and m c h = 0.35, for which R^ = 3.76 and -RJ3k = 9-5. Because of the effects of 
systematic errors in the determination of both .Rmk and i?bdj it is not possible to give 1 sigma error 
estimates for 7 and m c h. If we assume that either -RJ^K or ^bd nas * ne ^ es * estimate of its value 
(9.5 and 3.76, respectively) while the other variable has its maximum error, then 7 lies in the range 
0.42-0.58 and m c h lies in the range 0.30-0.43. As stated above (see eq. 24), the turnover mass m pea k 
of the STPL function is in general less than m cn . Figure S] also shows that m pea k = 0.27 ± 0.04. 
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3 3.5 4 4.5 5 5.5 6 6.5 

Fig. 3. — Position of the observed sample in the list of the 1000 synthetic samples when arranged 
in decreasing order of KS-test significance levels. The numerical experiment was repeated for R at 
0.0025 intervals. For the KS-tests a Gaussian cumulative distribution with a = 0.3 was used. 
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Fig. 4. — The isocurves Rmk = 8, 9.5 and 11 (dashed); the isocurves Eg* = 3.37, 3.76 and 4.04 
(continuous); and the isocurves m pea k = 0.23, 0.27 and 0.31 (grey) for the STPL function with 
T = 1.35. The circle indicates the standard STPL individual-star IMF V'stpl, which has parameter 
values (r = 1.35, 7 = 0.51 and m ch = 0.35) that fulfill the constraints = 3.76 and i? M K = 9.5. 
The errors in the results derived for ^stpl m Table 2 are the maximum errors when either -RjJJk 
or R™£ has the best estimate of its value (9.5 and 3.76, respectively). 
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4. Comparison and Characterization of IMFs 

Table 2 and Figures [2k., b,c summarize the results for the various IMFs described in § §2.2l [231 
In the table, we separate the results for individual (top) and system or effective (bottom) IMFs. 
Since the Salpeter IMF does not include brown dwarfs, we set the lower limit m.£ on this IMF with 
Rmk- The results for the Kroupa, Chabrier, Padoan & Nordlund, Hennebelle Sz ChabrieiEl and De 
Marchi et al IMFs correspond to the parameter values suggested by the authors. For comparison 
purposes we estimate that the effective STPL IMF, "0stpL' * s characterized by parameters 7 = 0.57 
and m c h = 0.42. These parameters are obtained by assuming that the effective peak mass is a 
factor 1.25 larger than the peak mass for individual stars (as Chabrier found) and that = 4.6. 

Figure [2^t shows that the Chabrier individual-star IMF differs by less than 30% from the 
STPL individual-star IMF for m > 0.03. It is is not surprising that both IMFs agree for m > 0.1 
since both rely on the RGH02 data, although the STPL IMF also makes use of data on brown 
dwarfs that was not available to Chabrier. If brown dwarfs are omitted from the IMF, then the 
agreement between the IMFs is even better (jl5%). Although the Chabrier (2003b) IMF ipco3 has 
been updated by Chabrier (2005), we include it in Figured since recent works continue to use it. 
The Kroupa individual-star IMF differs from V'stpl by l ess than a factor 1.6 for m < 1.4, but since 
it has r = 1.7 the difference at large masses increases substantially. Figure [2b shows the effective 
system IMFs. The STPL form of this IMF is based on Chabrier's (2005) result that the peak mass 
of the effective IMF is 1.25 times that for the individual-star IMF, so it is natural that the two 
are in good agreement, except again for m < 0.03. The Kroupa effective IMF differs by less than 
a factor 1.8 with respect to V'stpl ™ the entire mass range, but is much larger than all the other 
IMFs at low masses. The MPPZ IMF differs significantly at high masses, but it must be borne in 
mind that the high-mass portion of this IMF is based on data that extend only to 10 Mq ; in the 
mass range 0.04 < m < 10, the difference with V'stpl ^ s ^ ess than a factor 1.55. The Scalo (1986) 
IMF differs from V'stpl by less than a factor 2 in the mass range 0.1 - 3.4; note that the difference 
for m > 3 would be reduced if it were boosted by l/b(to) (see Paper II) and the factors mentioned 
in E j3.ll (Parravano et al. 2003, 2009). Figure [2J; shows that the theoretical Padoan-Nordlund IMF 
agrees well with the STPL IMF (and the Chabrier IMF) for m > 0.05. The Hennebelle-Chabrier 
IMF differs by up to a factor 1.8 from the STPL IMF in this mass range, but the data comparisons 
we have made in this paper do not allow us to determine which is more accurate. Note that HC 
IMF provides an example of an IMF that is flatter at intermediate masses than at high masses. It 



2 The Padoan-Nordlund and the Hennebelle-Chabrier theories provide the core mass function (CMF) — i.e., the 
mass function of gravitationally bound regions in a molecular cloud. Observations show that the CMF is similar 
to the IMF, but shifted up by a factor of about 3 (McKee & Ostriker 2007). Hennebelle & Chabrier (2008, 2009) 
therefore assumed that each core is three times the mass of the star(s) that form within it. Padoan et al. (2007) took 
advantage of the scale-free nature of self-gravitating simulations to adjust their core masses directly to the stellar 
masses. In comparing these IMFs with the observed one, it is important to note that the observed system IMF is 
an effective system IMF, as discussed in the Introduction; the difference between the actual system IMF and the 
effective one is small compared to the uncertainties in the relation between the CMF and IMF, however. 
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must also be noted that the theoretical IMFs are for single values of the parameters, whereas the 
actual IMF is the result of averaging over a range of physical conditions. 

For each IMF, Table 2 first gives the M to K ratio, -Rmk> an d the brown dwarf ratio, R^d- The 
observed values were used to determine the parameters of the STPL IMF, so it is consistent with 
these observations by construction. For the individual-star IMFs, the Salpeter IMF also agrees 
with the observed .Rmk by construction, since its value of rri£ was chosen to do this; however, it 
includes no brown dwarfs. The Chabrier (2005) IMF is consistent with both the observed values of 
both -Rmk and R^d] the latter agreement is significant, since -Rbd is based on data published after 
Chabrier (2005). We did not include the Chabrier (2003b) IMF in the Table, but we note that it 
has Rmk = H-9 and i?bd = 1-92, which are far from the observed values. The Kroupa IMF has 
a value of -Rmk that is slightly too high and a value of -Rbd that is significantly too low. For the 
effective system IMFs, the Chabrier, PN and HC IMFs are all within the errors of the observed 
Rmk, but only the Chabrier IMF is within the errors for the brown dwarf fraction. Note that the 
MPPZ IMF, which also has the STPL form, has a value of -Rmk that is far from the observed value. 

Next, Table 2 gives the mean mass of all stars, (m), the mean mass of main-sequence stars 
(i.e., omitting the brown dwarfs), (m) ms , and the peak mass, m pea k. As shown in Table 2, the 
peak masses for the PN, HC and Chabrier IMFs are in excellent agreement. The peak mass for 
the STPL IMF is somewhat higher than that for the Chabrier IMF, and that for the MPPZ IMF 
is somewhat lower. The Kroupa IMF has the lowest peak mass. According to the Chabrier (2005) 
IMF, the average observed mass of an individual main sequence star is (m ms ) = 0.69M Q ; the 
Salpeter IMF (with m e = 0.21) and the STPL IMF are consistent with this, but the Kroupa IMF 
gives a significantly smaller value, (m ms ) = 0.43. In fact, most of the individual-star quantities 
for the Kroupa IMF in Table 2 disagree with those for all the other IMFs. We interpret this as 
being due in part to his adoption of a steep slope, T = 1.7, for the high mass (m > 1) portion of 
the individual-star IMF. As noted above, Maiz Apellaniz (2008) has shown that the correction for 
binarity has only a small effect on the slope. 

The mass regimes in which disagreements are most significant are at the extremes, as might be 
expected. At the low-mass end of the IMF, Table 2 gives the fraction of stars that are born as brown 
dwarfs, Fbd = J\*{< m^/^T, where we adopt = 0.075. (Recall that we use rribd = 0.08 
only for comparison with the observationally determined -Rbd-) The fraction of stars that are born 
on the main sequence is F ms = 1 — Fbd- The brown dwarf number fraction depends slightly on rri£, 
which we assume to be 0.004. Increasing ni£ from 0.004 to 0.006 decreases the fraction of brown 
dwarfs by about 5% (and the mass fraction by only 1%). 

With the exception of the Kroupa, MPPZ and STPL IMFs, all the IMFs rapidly fall off at very 
low mass, so there is a clear distinction in the predictions of the number of low-mass brown dwarfs. 
The Chabrier (2005) IMF gives -Fbd = 0.22 for individual stars, whereas the STPL form of the IMF 
implies Fbd = 0.31. (We note parenthetically that the Chabrier 2003b IMF for individual stars has 
Fbd = 0.47, which is even larger than the Kroupa value Fbd = 0.39; however, as pointed out above, 
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the Chabrier 2003b IMF is inconsistent with observation.) The Chabrier (2005) and STPL IMFs 
differ significantly only for m < 0.03, where as yet there is little data. For the system IMFs, the 
Chabrier IMF implies Fbd = 0.17. The MPPZ IMF is in good agreement with this, but the PN 
and HC IMFs are below the Chabrier value by a factor of about 1.6. The STPL IMF has a value 
for the system brown dwarf fraction that is a factor of about 1.5 times greater than the Chabrier 
value; the Kroupa effective IMF has a brown dwarf fraction more than twice the Chabrier value. 
We note that brown dwarfs have a negligible fraction of the mass; for example, for ~ 1/3, the 
fraction of the mass in brown dwarfs is only ~ 2%. 
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Table 2. Characterization of IMFs 
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0.56 
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0.65 


0.77 


0.25 


0.17 
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85 
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14.4 
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0.15 
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0.57 
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a The IMF parameters are: V'SaiCT = 1.35, mi = 0.21, C = 0.167) for the Salpeter IMF; see 
eq. CL3]) with 73 = -1.7 and C = 0.076 for the Kroupa IMF; ^ c (m ch = 0.2, a = 0.55, T = 1.35, 
C = 0.315) in eqs. <^MB for the Chabrier IMF; VsTPL.(r = 1.35, 7 = 0.51, m ch = 0.35, C = 0.108) 
for the STPL with the constraints R™£ = 3.76 and -R^K = 9-5, the quoted errors give the values 
when -RjJJk: or ^bd nas ^ ne best estimate while the other variable has its maximum error. 

b System IMF: V>PN(r = 1.4, M = 10, m ch = 1, C = 0.101) for the Padoan and Nordlund IMF (i.e. 
a = 1.8); ^hc(?? = 0.4, M = 12, M* = y/2, m ch = 1.56, C = 0.218) for the Hennebelle and Chabrier 
IMF (i.e. a = 1.9). Effective IMFs: ^S86 from Scalo (1986) Table 7 for t = 12 Gyr and b(t Q ) = 1. 
For the effective Kroupa IMF see eq. (JEU) with 73 = -1.3 and C = 0.075; ^c(^ch = 0.25, a = 0.55, 
r = 1.35, C = 0.318— see eqs. HUES) for effective Chabrier IMF; Vmppz(P = 1.0,7 = 1.5, mch = 
0.15, C = 0.101) for the De Marchi et al. IMF; VstplO" = 1.35,7 = 0.57, m ch = 0.42, C = 0.146) 
for the STPL effective IMF with m poak = 0.34 and Rf d = 4.6. 

c Average mass of stellar objects (main-sequence stars and brown dwarfs) in the IMF. 

d Average mass of main-sequence stars in the IMF. 

e IMF turnover mass. 

f Fraction of stellar objects that are brown dwarfs (0.004 < m < 0.075). 
g Fraction of main-sequence stars with masses above = 8, normalized to 10~ 3 . 
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h Mass of stellar objects formed per high mass star. 
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For the high-mass end of the IMF, Table 2 gives two parameters describing stars that will 
explode as core-collapse supernovae, which have masses m > = 8. In general, we define 
Jx*h = A/"*(m > mh), the birthrate of massive stars, and Fh = AC/i/A^T; the fraction of all stars 
that are born with high masses. Table 2 gives the fraction of main-sequence stars that are born 
with high mass, F^ ms = Fh/F ms , and the mass ///, (in solar masses) of stars formed per high-mass 
star, 

»h = Ty- - (34) 
N*h b h 

where M* is the total star formation rate. The results in Table 2 show that the fraction of individual 
main-sequence stars that have high masses (m > 8) is F^ms — (7 — 8) x 1CP 3 , and the mass of stars 
formed per high-mass star is fih — (90 — 100)M Q . (We do not include the Kroupa individual-star 
IMF in this comparison because of the high value of T, as discussed above.) The fraction of the 
mass in high-mass stars is about 20% in each case. For the system IMFs, the Chabrier, PN and 
STPL forms all give i*X )ms ~ 9 x 10~ 3 and \ih ~ 90M Q . The value of \ih is smaller for the effective 
IMF than for the individual-star IMF because there are more stars with an effective mass above 
8 Mq than there are individual stars with such masses; stars with effective masses above 8 M & but 
actual masses less than that do not explode as supernovae. On the other hand, it is possible for 
one object in the effective IMF to produce more than one SN, an effect we did not allow for since 
it depends on the details of binary evolution. In any case, the individual-star value of fj,h should be 
used to estimate the rate of core collapse supernovae. (At present, however, this is just a matter 
of principle, since the uncertainties in each value are larger than their difference.) The Kroupa 
effective system IMF has a similar value of Hh but a somewhat lower value of the number fraction, 
Fh,ms- The HC IMF gives a somewhat higher fraction of high-mass stars and a correspondingly 
lower value of the mass formed per high-mass star. The MPPZ IMF has a value of the high-mass 
fraction that is about 1.5 times greater than the Chabrier value because of its relatively flat high- 
mass slope, r ~ 1; it should be borne in mind that De Marchi et al determined their mass function 
only for masses up to 1OM . 

In order to determine the effect of uncertainties in the high-mass slope, T, and the upper mass 
cutoff, m u , we have calculated the various quantities characterizing the IMF for the STPL IMF, 
holding mg and the low-mass constraints (Rmk = 9.5, i?bd = 3.76) fixed, but varying the high-mass 



Table 3. Characterization of V'STPl: dependence on T 



T 7 m ch m u (m) (m) ms m pcak F bd F^ ms /10 3 n h 

1.20 0.545 0.315 120+^ 0.55±g;g^ 0.77±g;g^ 0.26 0.30 10.6±g;| 74 

1.35 0.510 0.350 120^ 0.46lg;gJ 0.65±g;g^ 0.27 0.31 6.8±g;^ 97 

1.50 0.490 0.380 1201^ 0.40±g;gS 0.56±g;g? 0.28 0.31 4.4l°;° 131 
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slope (r = 1.2, 1.35 and 1.5) and the upper mass cutoff (m u = 80, 120 and 150). The upper mass 
cutoff has little effect on the quantities shown in Table 3. In particular, Hh = (m)/Ff l changes 
less than 1% with m u because both (m) and Fh increase with m u . The effects of varying V are 
shown in Table 3. The high-mass star fraction increases by somewhat more than a factor 2 as T 
decreases from 1.5 to 1.2; correspondingly, the mass of stars formed per high-mass star decreases 
by somewhat less than a factor 2. 

Star formation rates in galaxies are inferred from observations of high-mass stars. Since the 
shape of the IMF determines the ratio of the total mass of stars formed to the number of high- 
mass stars formed — in other words, \ih — the shape is essential for inferring the total star formation 
rate. It is striking that all the IMFs we have considered (with the exception of the individual-star 
Kroupa IMF, which has an anomalously steep high-mass slope) have a significantly greater fraction 
of high-mass stars, and as a result a significantly lower value of /ih, than does the Scalo (1986) 
IMF: for an 12 Gyr disk age, the Scalo IMF gives fih — 290, about 3 times greater than currently 
estimated. Correspondingly, as will be discussed in greater detail in Paper II, the current estimate 
for the Galactic star formation rate of M* = 1.3 ±0.2 Mq yr _1 (Murray & Rahman 2010) is several 
times less than earlier ones (e.g., 4 Mq yr _1 , McKee & Williams 1997), which were based on the 
Scalo IMF or variations of it. 



5. Observational Tests of the Shape of the IMF 

5.1. Comparison with Field Star Mass Function at Low Masses 

The IMF for individual stars, V'stplj can be directly compared with observed stellar mass 
functions only for samples in which binary systems are resolved or corrections are made for un- 
resolved systems. The STPL midplane mass distribution (continuous curve in Fig. [5^) is derived 
for a constant SFR during the last 11 Gyr and the scale height 2H{m) from Scalo (1986). The 
mean SFR during the disk age is adjusted to match a surface mass density of M dwarfs equal to 
11 Mq pc -2 (the Flynn et al. (2006) estimate for the thin disk, see Paper II). As stated in ^3.4[ 
most binary systems are resolved in the RGH02 stellar mass function for the local sample of stars. 
The histogram in Figure [5^ shows the RGH02 midplane mass distributions (see their figure 13a, b), 
one using the KTG m(My) relation (continuous histogram) and the other using the RGH m(My) 
relation (dashed histogram). These mass functions have upticks both at m ~ 0.1 and at m ~ 1.0. 
We assume that these apparent features are due to systematic effects such as uncertainties in the 
m{My) relations and the dependence of these relations on metallicity. From integration of the 
RGH02-KTG mass distribution in Fig. [5^ and under the assumption that the scale height for M 
stars is 325 pc (Scalo 1986), we obtain a surface mass density of M dwarfs S m s(< 0.63) = 10.5 
Mq pc~ 2 . This value agrees with the ~ 11M pc -2 estimated by Flynn et al. (2006) for the 
surface mass density of M dwarfs in the thin disk. 

The effective IMF , V'stpL' * s compared in Fig. [5b with mass functions of distant field un- 
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resolved systems. Observations with the Hubble Space Telescope (HST) can be used to infer the 
PDMF of distant field stars. Zheng et al. (2001, hereafter ZFGBS01) studied main sequence stars 
with My > 8. There is no Malmquist bias in this study since HST can see all the Galactic M stars 
in this direction. In addition, since M stars have lifetimes substantially in excess of the Hubble 
time, the time-averaged IMF and the PDMF (per unit area of disk) are proportional for these 
stars. ZFGBS01 considered two color-magnitude relations to infer the stellar masses. Assuming 
solar metallicity (method CMR1), they inferred S ms (< 0.63) = 14.3 db 1.3 M & pc -2 , which is 
in agreement with previous estimates by Gould, Bahcall and Flynn (1996, 1997) using the same 
color-magnitude relation. On the other hand, assuming that the color-magnitude relation is a func- 
tion of height above the Galactic plane due to the vertical gradient in metallicity expected in the 
population of M stars (method CMR2), they inferred S ms (< 0.58) = 12.2 ± 1.6 M pc -2 . Note 
that the ZFGBS01 sample includes thick disk stars also and therefore their surface mass density 
of M dwarfs is expected to exceed the 11 M & pc~ 2 value estimated by Flynn et al. (2006) for the 
thin disk. We also note that ZFGBS01 presented the results for the CMR1 method in order to 
check the consistency with their previous work, but they regarded the metallicity-gradient model 
as more physically motivated and so more reliable. However, the significant difference between 
the PDMFs derived using the CMR1 and CMR2 approximations illustrates the uncertainties that 
can be introduced when the distances to objects are not available and a color-magnitude relation 
must be used. Regarding the effect of binaries, Gould et al. (1996) showed that the correction for 
binaries has only a small effect on the inferred surface density S^t- The open circles connected 
by dashed lines in Figure [5b show the midplane mass distribution from the ZFGBS01 PDMF with 
the CMR2 method and a scale height H = 325 pc. Parravano et al. (2006) fitted the STPL IMF 
to the CMR2 data and obtained an effective system IMF characterized by 7 = 0.8 and m c h = 0.24; 
by comparison, the parameters of the STPL effective IMF found here (which focuses on Rmk and 
Rbd) are 7 = 0.57 and m c h = 0.42. The peak mass for the Parravano et al. (2006) IMF is 0.225 Mq, 
whereas it is 0.34 Mq for the one here. 

Bochanski et al. (2010) considered a sample of ~ 15 million SDSS point sources spread over 
8,400 square degrees and obtained an effective mass function (continuous histogram in Fig. [5b) 
in the range 0.1 < m < 0.8 that is well described by a log-normal distribution with m pea k = 0.25 
and a = 0.28. This estimate of m pea k is in agreement with Chabrier (2005) but, as mentioned in 
§3.41 a is too small and the ratio R^[k = 13.8 is too high. The surface mass density of M dwarfs 
inferred from this mass distribution and a scale height H = 325 pc is E ms (< 0.63) = 8.1 Mq pc~ 2 . 
Bochanski et al. (2010) have also estimated the mass function of single stars (not shown in Fig. (5b) 
by correcting for unresolved binaries, but their color-magnitude relation has not been corrected for 
possible metallicity changes with position. This individual-star mass function can be described by 
a log-normal distribution with m pe ak = 0.18 and a = 0.34, and the associated surface mass density 
of M dwarfs is S ms (< 0.63) = 10.0 Mq pc~ 2 . In view of the uncertainties in the scale height H, 
the agreement with the Flynn et al. (2006) surface density (11M pc~ 2 ) is satisfactory. These 
mass distributions can be fitted with a STPL function to infer an individual-star or an effective 
IMF, as was done in Parravano et al. (2006) for the Zheng et al. (2001, CMR2) data. A STPL 
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with 7 = 0.89 and m c h = 0.24 {m pea k = 0.23) fits the constraint R%k = 13.8 from the effective 
mass function in Bochanski et al. (2010) and RtJ = 4.6 from §3.5. However, in addition to the 
limitations mentioned in §3.4|, Bochanski et al. (2010) consider that the systematic uncertainty in 
the luminosity function is dominated by the assumed CMR. Additionally we note that Bochanski et 
al. (2010) considered a CMR that depends on metallicity, but no metallicity correction was applied 
to the mass-Mj relation; their sample contains stars up to a height of 2 kpc, and such stars have 
lower metallicity than those in the solar neighborhood (Ivezic et al. 2008). Figure [Sb shows that 
the midplane effective mass function derived by Bochanski et al. (2010) and the effective STPL 
derived here agree within the errors for m < 0.4, but as the mass increases the Bochanski mass 
function becomes progresively lower; their mass function is Salpeter-like in the mass range 0.32-0.8. 



5.2. Shape of the Brown Dwarf PDMF 

The ideal test to the proposed IMF in the brown dwarf regime is to compare it with the 
mass distribution of field brown dwarfs, or with the average shape of the mass functions in a 
representative sample of clusters. Unfortunately, as explained in §3.51 at present both approaches 
have severe difficulties. Nevertheless, Allen et al. (2005) estimate that at a 60% confidence level 
the mass function for objects with masses between 0.04 and 0.10 Mq is dN*/dm oc m -(°-3±0-6) ^ 
corresponding to 7 = 0.7 ± 0.6. This estimate of 7 for nearby field brown dwarfs is consistent with 
our result for the effective system IMF, 7 = 0.57. Caballero et al. (2007) studied the young open 
cluster a Orionis and derived its substellar mass function down to 0.006 Mq . The masses were 
derived using the theoretical mass-luminosity relation of Baraffe et al. (2003). The main result 
from this study is that the mass function is continuous from low- mass stars to substellar objects, 
even below the deuterium-burning mass limit, with a power-law slope 7 = 0.4 ± 0.2 in the mass 
range 0.006 < m < 0.11, in agreement with our standard slope values. 

Current data are only now becoming adequate to measure the shape of the IMF in the brown- 
dwarf mass range, and with improving accuracy this measurement will provide a powerful discrim- 
inant among different forms for the IMF. For IMFs that are power laws at low masses, such as the 
Kroupa IMF, the MPPZ IMF and the STPL IMF, the fraction of brown dwarfs below a mass m is 
(eq. [25D 

Jbd(< m) _ rrf-ml 
F hd ml d - m < 

For example, the fraction of brown dwarfs with m < 0.03 Mq (below the mass range used to 
calculate i?bd) is 46%, 25% and 52% of the total number of brown dwarfs for the Kroupa, MMPZ 
and STPL individual-star IMFs, respectively (we used me = 0.004). IMFs with a log normal 
behavior at low masses predict smaller values for this ratio than either the Kroupa or STPL IMFs; 
for example, the individual-star Chabrier (2005) IMF has 30% of the brown dwarfs below 0.03 Mq. 
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log(m) log(m) 



Fig. 5. — Comparison of the midplane mass function corresponding to the STPL PDMF for a 11 
Gyr old disk with 11 M pc~ 2 in M dwarfs (see text) and the midplane mass distributions derived 
from field stars observations. The error bars beside the STPL label indicates the error from 2H, 
assumed to be 25%, and the uncertainty in the surface density of M dwarfs, taken to be in the range 
10 to 12 Mq pc~ 2 . The numbers beside the right hand ordinate gives the mass distribution per 
ln(m) units, a) The two individual-star mass distributions derived by RGH02 derived from nearby 
field stars (see text) and the distribution corresponding to the individual-star IMF V'STPL, which 
has parameter values (r = 1.35, 7 = 0.51 and m c h = 0.35). b) The midplane mass distributions 
corresponding to the effective IMF V'STPL) which has parameter values (r = 1.35, 7 = 0.57 and 
m ch = 0.42), and the effective mass distributions derived from distant field stars by Zheng et al. 
(2001) and by Bochanski et al. (2010). 
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6. Conclusions 

The stellar initial mass function (IMF) is fundamental to many branches of astronomy. Ob- 
servations suggest that when averaged over long times and large areas of the Galactic disk, it is 
universal in shape (Bastian et al. 2010). We argue that it should also be simple in shape, since there 
are no known processes that could imprint features on the IMF given the wide range of conditions 
under which stars form. We propose that the IMF can be represented by a smoothed two-power 
law (STPL), 

V'STPlim.) = Cm~ T {l - exp [-(m/m ch ) 7+r ] } (m e < m < m u ). 

This IMF has five parameters: The lower-mass cutoff, mi, which we set equal to the minimum 
fragment mass expected on the basis of opacity effects, 0.004 M & (Low & Lynden-Bell 1976; 
Whitworth et al. 2007); the upper mass cutoff, m u , which appears to be about 12OM0 for individual 
stars ( §3, 2ft : the high-mass slope parameter, T, which we take to be the Salpeter value of 1.35 ( §3. lft : 
the low-mass slope, 7; and the characteristic mass, m c h, which is comparable to the mass at which 
the IMF peaks, 

m pea.k- In order to determine these latter two quantities, we use integral properties 
of the IMF that are set by observation: (1) the ratio -Rmk of the number of stars (mainly M-dwarfs) 
in the range 0.1 — 0.6M Q to the number of stars (mainly K dwarfs) in the range 0.6 — 0.8Mq, as 
inferred from the Reid et al. (2002; RGH02) sample of nearby stars ( §3.4p : and (2) the ratio i?bd 
of the number of stars in the mass range 0.08 — 1M divided by the number of objects (mostly 
brown dwarfs) between 0.03 and O.O8M ( §3.5p . The RGH02 data are for individual stars, but the 
brown dwarf data must be corrected for unresolved binaries. With the observed values of Rmk and 
i?bd) we infer 7 = 0.51 and m c h = 0.35 for the STPL individual-star IMF; the corresponding peak 
mass is m pea k = 0.27. For these values of the parameters, the normalization factor is C = 0.108. 
The STPL form for the effective system IMF is found by using the value of i?bd found for systems 
in §3.51 and assuming that the peak mass for systems is a factor 1.25 larger than the peak mass 
for individual stars, as Chabrier (2005) found. The resulting parameter values are 7 = 0.57 and 
n^ch = 0.42; the corresponding peak mass and normalization factor are respectively m pea k = 0.34 
and C = 0.146. 

A STPL form for stellar mass functions was first proposed by Paresce & De Marchi (2000), 
who applied it to globular clusters. De Marchi & Paresce (2001) showed that this form applied to 
young Galactic clusters as well, and drew the important conclusion that the characteristic mass of 
Galactic clusters increases with the age of the cluster. Our approach in studying the mass function 
is complementary to theirs: They have focused on analyzing individual clusters, whereas we have 
used data on nearby stars to infer the shape of the IMF for low-mass main sequence stars and a 
broad range of data to infer the high-mass slope of the IMF. We have also considered a broader 
mass range, using data from 0.03 Mq to over 100 Mq, whereas they have focused on the mass 
range between 0.1 Mq and 10 Mq. For their average effective mass function, they find T ~ 1, 
7 ~ 1.5, and m c h = 0.15; the resulting values for -Rmk (14.4) and i?bd (5-7) are not consistent with 
the values we have determined. The reasons for this discrepancy are not entirely clear. We note 
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that the parameter 7 determines the number of substellar objects; we use observations of the ratio 
of the number of substellar objects to the number of low-mass stars (i?bd) to infer 7, whereas De 
Marchi et al. (2010) do not include objects less than 0.1 in their analysis. 

The STPL IMF proposed here agrees well with the Chabrier (2005) IMF (i.e. to within 30% 
for m > 0.03). Both are based on the RGH02 data, but the STPL IMF also makes use of data 
on brown dwarfs that was not available in 2005. As shown in Table 2, the mean stellar mass is 
(m) ~ 0.5; if only main sequence stars are included, the mean mass is (m) ms ~ 0.7. This agreement 
applies to both the individual-star and effective system IMFs; the latter is not surprising, since we 
used Chabrier's (2005) result on the ratio of the peak masses in converting the individual-star IMF 
to the effective system IMF. Chabrier's earlier version of the individual-star IMF (Chabrier 2003b) 
is inconsistent with both the value of -Rmk determined from local stars and with the value of i?bd 
determined from very young clusters. 

As shown in figured, the Salpeter form of the individual-star IMF, which had its lower mass 
limit m,£ adjusted to fit the observed value of Rmk, is above the STPL IMF at all masses since 
its normalization factor (C = 0.167) is 1.55 times the corresponding STPL value. The Salpeter 
and STPL IMFs agree within 30% for m > 0.4 when the Salpeter IMF is reduced by a factor 1.55 
or when the IMFs are normalized to produce a unit mass. At large masses there are significant 
disagreements between the Kroupa individual-star IMF and the others since Kroupa's high mass 
slope (r = 1.7) is substantially larger than others, but for m < 1.4 it agrees to within a factor 1.6 
with the STPL IMF. 

Figure [2)3 shows the theoretical system IMFs. Both the Padoan & Nordlund (2002) and 
Hennebelle & Chabrier (2008, 2009) forms for the IMF are in reasonable agreement (within 30%) 
with the STPL and Chabrier IMFs for m > 0.05, except for rn > 1 where Hennebelle & Chabrier 
IMF exceed STPL IMF up to a factor 1.8. They are both within the errors of the measured value 
of -Rmk> but they underpredict the number of brown dwarfs (i?bd is too high). The success of these 
theoretical forms for the IMF in matching the data is impressive, particularly in view of the fact 
that the actual IMF is the result of stars forming in a wide range of physical conditions, whereas 
the theoretical values are for only one set of conditions. 

The mass range in which the disagreement among the IMFs is greatest is the brown dwarf 
regime, which is also the mass range in which the data are most limited. The STPL, MPPZ and 
Kroupa forms for the effective IMF assume that the very low mass range is described by a power 
law, whereas Chabrier assumes a log normal for that mass range. Both theoretical IMFs (PN and 
HC) have a log-normal form built in from the log normal form for the density spectrum found in 
numerical simulations. A sensitive measure of the IMF in the brown dwarf regime is the fraction 
of brown dwarfs predicted to be below some mass; for example, under the assumption that brown 
dwarfs extend down to 0.004 Mq, the STPL and Kroupa forms for the IMF predict that about 
half of all brown dwarfs have masses below 0.03 Mq, whereas the Chabrier form predicts that only 
about 27% of brown dwarfs have such a low mass; PN and HC forms predict even lower values, 
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14% and 15% respectively. Since the formation of brown dwarfs requires unusual conditions, the 
IMF of brown dwarfs is a powerful test of star formation theories (McKee & Ostriker 2007). 

The IMFs considered here agree reasonably well on the proportion of high-mass stars (m > 
m h = 8), which will produce core-collapse supernovae. For the individual-star IMF, the mass of 
stars formed per high-mass star is fih — 90 — 100. (The effective system IMF has a somewhat lower 
value of Hh because it includes stars of mass m < 8 in binaries with effective masses above 8 M Q .) 
This value of fih is significantly smaller than that for the Scalo (1986) IMF, which has fih = 280. 
As will be discussed in greater detail in Paper II, the result is that current estimates of the star 
formation rate in the Galaxy are several times lower than those in previous decades. 

We acknowledge support from the NASA Astrophysical Theory Program in RTOP 344-04-10- 
02, which funded The Center for Star Formation Studies, a consortium of researchers at NASA 
Ames, University of California at Berkeley, and University of California at Santa Cruz. AP was 
supported as a Senior Associate for part of this research by the National Research Council and 
by the Universidad de Los Andes by CDCHT project C-1653-09-05-B. The research of CFM was 
supported in part by NSF grant AST-0098365 and by the Groupement d'Interet Scientifique (GIS) 
"Physique des deux infinis (P2I)." 

A. Binary Correction 

Recall from §1 that there are three forms of the IMF, and correspondingly three ways to 
count stellar and substellar systems: Af* nd is the number of individual stars, counting two stars for 
each binary, three for each triple, etc; M sys is the number of stellar and substellar systems in which 
binaries, etc., are counted as single objects and characterized by their total mass (or primary mass); 
and Af eS is the number of stellar and substellar systems in which each system is characterized by 
an effective mass, m e g, determined from unresolved observations in which multiple systems are 
assumed to be single. In the text, we use two observed ratios to fix the form of the IMF: Rmk, 
which is roughly the ratio of M stars to K stars, and R^, which is the ratio of the number of low- 
mass stars to readily observable brown dwarfs (actually, stars in the mass range 0.03 — 0.08 Mq). 
Current data allow one to determine Rmk for individual stars, but the data for i?bd do not resolve 
binaries; the data therefore determine i^d- To relate the observed R^ to the required R^, we 
first need to account for binaries and then allow for the difference between the effective mass and 
the actual mass. 

We denote the mass range for low-mass stars (0.08 < m < 1) by the subscript and the mass 
range for readily observable brown dwarfs (0.03 < m < 0.08) by the subscript "bd." For example, 
Af^f is the number of low-mass stellar systems with effective masses in the low-mass range and 
A/^f is the number of substellar systems with effective masses in the range 0.03 — O.O8M . The 
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observed ratio i?bd is based on infrared observations that did not resolve binaries 



-Rbd — R\ 



,eff _ 
•bd = 



Aff(0.08 < m cS < 1) _Mtf 



(Al) 



AC cff (0.03 < m cS < 0.08) A/?f 



whereas we need 



-"1 



>ind 

bd = 




(A2) 



to fix the form of the IMF. 



A.l. Ratio of Low-Mass Stars to Low-Mass Systems, M'^ d /M\ 



-sys 



First we consider the ratio of the number of individual low-mass stars to the number of low-mass 
systems. We characterize binaries by the mass of the primary, m\. The subscript "6" denotes binary 
stellar systems; thus, Mb,t is the number of binary stellar systems in which both the primary and 
the secondary stars are in the low-mass range. The number of stellar systems in which the mass of 
the star, for single systems, or the mass of both the primary and the secondary, for binary systems, 
is in the low-mass range is denoted MJf 5 - (Thus, systems are placed in the low-mass category 
based on the individual stellar masses, not on the total mass of the system.) The binary fraction 
for low-mass stellar systems is then /& £ = A/& e/J^J ■ Note that brown-dwarf (BD) companions 
of low-mass stars are not counted in fa £\ these will be discussed below. For the moment, let us 
neglect the number of low-mass companions of stars with masses exceeding 1 M & (which we shall 
term "supersolar stars"), we have 



According to Lada (2006) the PDMF of locally observed stars has a single star fraction of 2/3. 
Since most of the stars in the PDMF are less than 1 M Q , this corresponds to a single star fraction 
below 1 Mq of 2/3, or a binary fraction £ ~ 1/3. 

We now determine the correction for low-mass companions of supersolar stars. For a flat 
distribution of secondary to primary mass (Mazeh et al. 2003), one can show that about half 
the supersolar binaries have a low-mass companion. We obtain an upper limit on the effect of 
these low-mass companions by assuming a binary fraction of unity for supersolar stars; then the 
number of low-mass companions of supersolar stars is 1 /2 the number of supersolar systems, and 
the number of supersolar stars is 3/2 times the number of supersolar systems. The ratio of low-mass 
companions of supersolar stars (denoted by the subscript "ss") to low-mass stars is then about 



where we obtained the ratio of the number of supersolar stars to low-mass stars in the individual- 
star IMF by iteration. Thus, the number of low-mass stars in equation (|A3|) should be increased 
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by about 5%. We conclude that 



Mlf * 1.05 x | ACf = 1.4QACT- (AS) 



A.2. Correction for the Effective Mass, JVT/A/] 



cff 



Next, we must estimate the ratio of the number of low-mass stellar systems, ACj , to the 
number of stellar systems with low effective masses, AfjfF • Fortunately, this correction too is very 
small. We approximate the luminosity in a photometric band as a power-law in mass, 



L 



pb 



X 



m 3 P b _ 



(A6) 



The results of Delfosse et al. (2000) show that j p b — 2 for 0.1 < m < 0.7 for both the J-band 
and the K-band (their more complicated fits are needed to extend the fits over a slightly larger 
mass range). For a low- mass binary with a mass ratio q = m^/mi, where is the mass of the 
secondary, the luminosity of the system corresponds to an effective mass 



m eS = mi(l + q 



2U/2 



(A7) 



We estimate the effect of the effective mass by assuming that all the binaries have the typical value 
of q. For low-mass systems, this is about | according to Mazeh et al. (2003), corresponding to 
mi = 0.89 m e g. The effective number of low-mass stellar systems, Af^f, omits the binaries in the 
mass range 0.89 — 1 Mq and includes the binaries in the mass range just below 0.08 Mq. We obtain 
an upper limit on the ratio N^ s /M^f by omitting the latter: 



ACf ^ KT - A4(o.89 < mi < i). 



(A8) 



The binary fraction for solar mass stars was inferred to be 0.57 by Duquennoy & Mayor (1991). 
(The Lada (2006) binary fraction used above for m < 1 is smaller because it includes less massive 
stars.) Using equation (|A3p . we then obtain 



■cff 



A/?f 



i 



i 



0.57AC ys (0.89 < mi < 1) 



0.57(1 + / M ) A/"i nd (0.89 < mi < 1) 



1.57 



A/l? d 



(A9) 
(A10) 



The determination of the ratio A/^ nd (0.89 < mi < 1)/Af^ d must be done iteratively; we find that 
it is about 0.03 so that the correction in going from the effective mass to the primary mass results 
m an increase of at most about 1%, implying A/Jf ~ AQf . We conclude that Af'^ d ~ 1-4A/J/ . 
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A.3. Effect of Brown Dwarfs in Binaries, A^/AT^ 

We turn now to the effect of brown dwarfs in binaries, which enters the denominator of i?bd- 
Very low- mass (VLM) binaries are defined as those with primary masses less than 0.1 M@. Reid et 
al. (2006) and Burgasser et al. (2007) find VLM binary fractions of 24^% and 22+^%, respectively. 
Since the binary fraction tends to decrease as the mass decreases (Burgasser et al. 2007), we adopt 
a binary fraction for BD-BD binaries of 20%, recognizing that this is uncertain. Since BD-BD 
binaries tend to have similar masses (Burgasser et al. 2007), we assume that all the BD secondaries 
are in the target mass range of 0.03 — 0.08 M@. 

BD secondaries of H-burning stars are comparable in number to the BD secondaries of BD 
primaries. Burgasser et al. (2007) summarize the data on BD companions of low-mass stars. 
Imaging surveys indicate that the number of BD secondaries is about 1/3 that of hydrogen-burning 
secondaries; radial velocity surveys, which are sensitive to smaller separations, give a somewhat 
smaller fraction. We adopt a value fb,e/3 for the fraction of low-mass systems with BD secondaries, 
recognizing that this too is uncertain. Then the total number of brown dwarfs is the number of 
systems with BD primaries, plus the number of BD companions of BD primaries, plus the number 
of BD companions of H-burning stars: 

Kf = AC +M>d,6 + \h,eKF- (All) 

Keep in mind that in order to conform to the convention adopted by Andersen et al. (2008), we 
are including stars with masses between the H-burning limit ~ 0.075 M & and 0.08 M Q in Abd- We 
have just seen that Abd,6 — 0-2A/"^ s , so we can solve this equation to find 



pind V _ *t I * _ *t J "bd ^ bd ( 



since the binary fraction of low-mass stars is fb,t — § (Lada 2006) and we previously found that 
the total number of low-mass stars is about 1.4 times the number of low-mass systems. 

In principle, there could be a correction in going from the system counts to the effective system 
counts (i.e., from A/^ s to A^f), but as in the case of the low-mass stars, we expect this correction 
to be small. The simple approach adopted for estimating this effect for low-mass stars (starting 
with eq. IA6|) does not work for brown dwarfs, since spectroscopic information is also used in 
estimating the mass of the brown dwarf; furthermore, for field brown dwarfs, the estimate of the 
mass is coupled to an estimate of the age. Since this correction is small, and we earlier found that 
A^ s ~ Af^f , we have ~ R\/d ■ According to the estimation in §3.51 R^ = 4.6lgg 5 . Equation 
(|Al2l> then implies R^ ~ 3.76±£f. 
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